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Abstract 

We investigate transverse hadron spectra from proton+proton, deuteron+Au and 
Au+Au collisions at yfs = 200 GeV and =17.3 GeV within the Hadron-String- 
Dynamics (HSD) approach which is based on quark, diquark, string and hadronic 



degrees of freedom as well as Pythia calculations for the high p± spectra. The 
comparison to experimental data on transverse mass spectra from p+p, d+Au and 
Au+Au reactions shows that pre-hadronic effects are responsible for both the hard- 
ening of the spectra for low transverse momenta (Cronin effect) as well as the sup- 
pression of high p± hadrons. The interactions of formed hadrons are found to be 
negligible in central Au+Au collisions at y/s = 200 GeV for p± > 6 GeV/c, but have 
some importance for the shape of the ratio Raa at lower p± values (< 6GeV/c). 
The large suppression seen experimentally is attributed to the inelastic interactions 
of 'leading' pre-hadrons with the dense environment, which could be partly of col- 
ored partonic nature in order to explain the large attenuation seen in most central 
Au+Au collisions. 
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1 Introduction 



The phase transition from partonic degrees of freedom (quarks and gluons) 
to interacting hadrons - as expected in the early universe shortly after the 
"big bang" - is a central topic of modern high-energy physics. In order to 
understand the dynamics and relevant scales of this transition laboratory ex- 
periments under controlled conditions are presently performed with ultra- 
relativistic nucleus-nucleus collisions. Hadronic spectra and relative hadron 
abundancies from these experiments reflect important aspects of the dynam- 
ics in the hot and dense zone formed in the early phase of the reaction. In fact, 
estimates based on the Bjorken formula [1] for the energy density achieved in 
central Au+Au collisions suggest that the critical energy density for the for- 
mation of a quark-gluon plasma (QGP) of 0.7 to 1 GeV/fm 3 [2] is by far 
exceeded in the initial phase for a couple of fm/c at Relativistic Heavy Ion 
Collider (RHIC) energies [3], but it is still the aim to unambiguously identify 
the formation and properties of this new phase. 

Presently, transverse momentum (or transverse mass) spectra of hadrons are 
in the center of interest. Here the suppression of high transverse momentum 
hadrons is investigated in Au+Au reactions relative to p+p collisions at RHIC 
energies of ^/s = 200 GeV [4] , since the propagation of a fast quark through a 
colored medium (QGP) is expected to be different from that in cold nuclear 
matter as well as in the QCD vacuum due to the energy loss of induced gluon 
radiation [5,6,7,8,9]. In fact, the PHENIX [10], STAR [11] and BRAHMS [12] 
collaborations have reported a large relative suppression of hadron spectra for 
transverse momenta above p± ~ 3 • • -4 GeV/c which might point toward the 
creation of a QGP, since this suppression is not observed in d+Au interac- 
tions at the same bombarding energy per nucleon [12,13,14]. Accordingly, the 
experimental observations are qualitatively in line with expectations from per- 
turbative Quantum Chromo Dynamics (pQCD) [15,16,17,18,19], but do not 
support the idea of initial state gluon saturation [20]. But, at present, it cannot 
be ruled out that a larger fraction of the suppression seen in central Au+Au 
collisions might be also due to hadronic final state interactions as suggested 
in Ref. [21]. It is the aim of this work to quantify the amount of hadronic final 
state interactions in d+Au and Au+Au collisions at RHIC energies. 

We employ the HSD transport model [22,23,24] for our study. This approach 
takes into account the formation and multiple rescattering of formed hadrons 
as well as unformed 'leading' pre-hadrons and thus superseeds the incoher- 
ent summation of individual p+p collisions. Such transport calculations allow 
to study systematically the change in the dynamics from elementary baryon- 
baryon or meson-baryon collisions to proton-nucleus reactions or from periph- 
eral to central nucleus-nucleus collisions in a unique way without changing any 
parameter. In the HSD approach nucleons, A's, N*(1440), N*(1535), A, E and 
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E* hyperons, S's, H*'s and JTs as well as their antiparticles are included on the 
baryonic side whereas the 0~ and 1~ octet states are included in the mesonic 
sector. Inelastic hadron-hadron collisions with energies above \/s ~ 2.6 GeV 
are described by the Fritiof model [25] (including Pythia v5.5 with Jet- 
set v7.3 for the production and fragmentation of jets [26]) whereas low energy 
hadron-hadron collisions are modeled in line with experimental cross sections. 
We mention that no explicit parton cascading is involved in our transport cal- 
culations contrary to e.g. the AMPT model [27]. 

A systematic analysis of HSD results and experimental data for central nucleus- 
nucleus collisions for 2 • • ■ 160 AGeV has shown that the spectra for the 'longi- 
tudinal' rapidity distribution of protons, pions, kaons, antikaons and hyperons 
are in reasonable agreement with available data. Only the pion rapidity spectra 
are slightly overestimated from AGS to SPS energies [28] which implies, that 
the maximum in the K + /n + ratio at 20 • • • 30 AGeV - seen in central Au+Au 
(Pb+Pb) collisions [29] - is missed. For a comparison of HSD calculations 
with experimental data at RHIC energies we refer the reader to Ref. [30]. 

In this work we concentrate on the transverse momentum dynamics and es- 
pecially on the very high momentum tail of the hadron spectra. In order to 
describe these high p± spectra, we use the Pythia v6.2 event generator [31] 
for nucleon-nucleon collisions, which describes the high transverse momentum 
spectra of peripheral nucleus-nucleus collisions from RHIC at y/s = 130 GeV 
very well [21]. We recall, that within this framework also experimental dilep- 
ton and direct photon spectra at SPS energies [32] or single lepton spectra at 
RHIC energies from open heavy quark production [33] are well reproduced. 



2 Transport calculations versus experimental data 

We start with hadron production in p+p collisions at the invariant energy 
y/s = 200 GeV. The comparison of the calculations for the transverse momen- 
tum spectra of charged hadrons from p+p collisions at midrapidity with the 
experimental data of the PHENIX and STAR collaborations [34,11] is shown 
in Fig. 1, which demonstrates that the description of hadron production in 
the elementary reaction is quite well under control (see also Ref. [21]). 

2. 1 Preliminaries 

Before coming to proton+nucleus or deuteron+nucleus collisions it is impor- 
tant to specify the concept of "leading" (pre-) hadrons and "non-leading" or 
"secondary" hadrons since this separation is of central importance for the re- 
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Fig. 1. The invariant cross sections for the production of neutral pions and charged 
hadrons in p+p collisions (\/s = 200 GeV) at midrapidity as a function of transverse 
momentum. The Pythia (version v6.2) calculations (solid lines) are compared to 
the experimental data from RHIC (open circles: PHENIX [34], filled squares: STAR 
[11]). The charged hadron spectra are divided by a factor of 100 for clarity. 

suits to be discussed below. To this aim we recall that in the string picture 
- in a high energy nucleon+nucleon collision - two (or more) color-neutral 
strings are assumed to be formed, which phenomenologically describe the low 
energy "coherent" gluon dynamics by means of a color electric field, which is 
stretched between the 'colored' ends of each single string (cf. Fig. 2a for an il- 
lustration). The latter string ends are defined by the space-time coordinates of 
the constituents, i.e. a diquark and quark for a "baryonic" string or quark and 
antiquark for a "mesonic" string. These constituent quarks, diquarks or anti- 
quarks are denoted as "leading" quarks. In line with Ref. [35] these constituent 
'partons' are assumed to pick up (almost instantly) an anti-colored parton 
from the vacuum and achieve color neutrality. These color-neutral objects - 
containing the string ends - are denoted as "leading pre-hadron" (cf. Fig. 2b 
for an illustration of the string fragmentation in space and time in the rest 
frame of the string.) 

In principle, a leading pre-meson or pre-baryon is able to hold up to 2 (3) 
leading quarks, but in the case of interest, i.e. midrapidity transverse momen- 
tum spectra in heavy-ion collision at \fs = 200 GeV, the majority of leading 
mesons/hadrons consists only of one leading quark and some secondary (di-) 
quarks, as will be discussed below. We note, that this changes dramatically 
for rapidity regions closer to those of the projectile or target or different y/s 
values, as discussed below. 

The time that is needed from the instant of the 'hard' pQCD collision t = 
to the formation of qq or qqqq pairs from the vacuum (ti > t ) and for the 
hadronization of the fragments (£2 = T/) we denote as formation time 77 in line 
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with the convention in transport models. For simplicity we assume (in HSD 
and also for the high p±_ partons), that the formation time is a constant ij 
in the rest frame of each hadron and that it does not depend on the particle 
species. Other estimates (cf. [21,36,37]) show a more complicated structure 
and particle dependence (especially for the light pion), but this issue will not 
be discussed in detail here; nevertheless, we will show results for different 
formation times and alternative scenarios below. In principle one expects also 
a distribution in the formation times, however, we here address only the mean 
value of such an 'unknown' distribution. We recall, that due to time dilatation 
the formation time tf in any reference frame is then proportional to the energy 
of the particle 

t f = 1 -T f = ^--Tf. (1) 



The size of 77 can roughly be estimated by the time that the constituents 
of the hadrons (with velocity c) need to travel a transversal distance of a 
typical hadronic radius (0.5 • • -0.8fm). Furthermore, since after the formation 
of a color neutral pre-hadron (t > ti) the Fock-components of this object are 
widespread in mass according to the uncertainty relation in energy and time, it 
will take a time of AE^ 1 to form a specific hadronic state of well defined radial 
excitation, where AE denotes the level distance of the hadronic spectrum. 

We assume in our transport simulations that hadrons, whose constituent 
quarks and antiquarks are created from the vacuum in the string fragmen- 
tation (at times ti), do not interact with the surrounding nuclear medium 
within their formation time (cf. Fig. 2c for an illustration). For the leading 
pre-hadrons, i.e. those involving (anti-)quarks or diquarks from the struck nu- 
cleons, we adopt a reduced effective cross section cri cad during the formation 
time Tf and the full hadronic cross section later on. Since this rough approxi- 
mation is subject of current debate [35] we will also present alternative models 
for the cross section of the leading pre-hadrons (see below). The 'default' con- 
cept is illustrated in more detail in Fig. 2c where a 'star' signals a possible 
collision while crossing lines imply that no collision might occur. 

Due to time dilatation light particles emerging from the middle of the string 
can escape the formed 'fireball' without further interaction, if they carry a 
high momentum relative to the rest frame of the fireball. We note that these 
arguments are strictly valid only in case of a constant formation time Tf in 
the rest frame as described above and implemented here. However, hadrons 
with transverse momenta larger than 6GeV/c predominantly stem from the 
string ends and therefore can interact directly with a reduced cross section 
(see below). 
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Fig. 2. Upper Part: Illustration of the formation of two strings in a N+N collision. 
The formed strings are color-neutral in total, however, carry (anti-)color at the 
string-ends, which here are displayed as diquark (d) - quark (q) pairs. The hatched 
areas denote nonperturbatively interacting partons. Middle Part: The space-time 
evolution of a string S in its rest frame. Whereas the 'leading' (di-)quarks (repre- 
sented by the outer lines) form 'leading' pre-hadrons almost instantly, secondary 
qq pairs are excited in the non-perturbative vacuum (after some time delay t±) and 
recombine to secondary hadrons (shadowed double lines). The time from the initial 
creation of the string to the final formation of secondary hadrons is denoted as Tf . 
The distribution in the formation time is presently unknown and for convenience 
characterized by a single formation time tj for all hadrons (dotted hyperbola). 
Lower Part: Whereas 'secondary hadrons' are not allowed to scatter during their 
formation time, they interact with the full hadronic cross section after being formed. 
The 'leading' color neutral pre-hadrons, however, carrying a constituent quark from 
the struck hadron, may scatter (in-)elastically with other leading pre-hadrons as 
well as formed secondary hadrons with a reduced cross section. In the present real- 
ization of HSD this reduced cross section is determined in line with the constituent 
quark model (see text). 
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2.2 Numerical implementation 



For the production and propagation of hadrons with high transverse mo- 
mentum (> 1.5GeV/c) we employ a perturbative scheme as also used in 
Refs. [38,39,40,41] for the charm and open charm degrees of freedom. Each 
high p± hadron is produced in the transport calculation with a weight Wi 
given by the ratio of the actual production cross section divided by the inelas- 
tic nucleon-nucleon cross section, e.g. 

x 

a NN W S ) 



In the transport simulation we follow the motion of the high p± hadrons within 
the full background of strings/hadrons by propagating them as free particles, 
i.e. neglecting in-medium potentials, but compute their collisional history with 
baryons and mesons or quarks and diquarks. For reactions with diquarks we 
use the corresponding reaction cross section with baryons multiplied by a fac- 
tor of 2/3. For collisions with quarks (antiquarks) we adopt half of the cross 
section for collisions with mesons and for the leading pre-hadron (formed) 
baryon collision a factor of 1/3 is assumed. The elastic and inelastic interac- 
tions with their fractional cross section are modeled in the HSD approach in 
the same way as for ordinary hadrons with the same quantum numbers via 
the Fritiof model [25] (including Pythia v5.5 with Jetset v7.3 for the 
production and fragmentation of jets [26]). 

The relative quark counting factors mentioned above might appear arbitrary 
and simplistic. However, this concept - oriented along the additive quark 
model - has been proven to work rather well for nucleus-nucleus collisions 
from AGS to RHIC energies [28,30,42] as well as in hadron formation and 
attenuation in deep inelastic lepton scattering off nuclei [43]. Especially the 
latter reactions are important to understand the attenuation of pre-hadrons or 
ordinary hadrons with high momentum in cold nuclear matter [35] . Our stud- 
ies in Ref. [43] have demonstrated that the dominant final state interactions 
(FSI) of the hadrons with maximum momentum, as measured by the HER- 
MES collaboration [44], are compatible to the concepts described above. This 
also holds for antiproton production and attenuation in proton-nucleus colli- 
sions at AGS energies [45]. Both independent studies point towards a hadron 
formation time Tf in the order of 0.4 • • • 0.8 fm/c. 

As mentioned above, a crucial question for the interaction dynamics is the frac- 
tion of leading pre-hadrons to secondary hadrons (or all produced hadrons) as 
a function of p± for the different hadron species. This information is directly 
extracted from the Pythia calculations for different hadrons by explictly trac- 
ing back the string fragmentation histories with all their details, i.e., begin- 
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ning with the incoming hadrons, following the hard interaction and linking 
the quark content of the final hadrons back to the partons at the initial string 
ends. The resulting fraction of leading pre-hadrons to all produced hadrons 
is displayed in Fig. 3 for N+N collisions at y/s = 200 GeV. We have checked, 
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Fig. 3. The ratio of 'leading particles' to 'all produced particles' (before particle de- 
cays) in N+N collisions (y/s = 200 GeV) at midrapidity for different particle classes 
as a function of transverse momentum within the Pythia description. (An isospin 
average in line with central Au+Au collisions has been applied.) For the region 
p± = 5 • ■ ■ 7GeV/c the theoretical error is of the order ±0.05 and to be neglected 
outside this region. 

that in Fig. 3 contributions with 2 (or more) leading particles are negligible 
and thus the above mentioned "leading" cross sections cri ea d as a fraction of 
the hadronic cross section of 1/2 for mesons and 1/3 for baryons are sustained 
by simple quark counting rules. 

In detail, one notices slight differences between pions and especially antibaryons 
(not shown), however, the fraction of leading pre-hadrons increases almost 
linearly with pj_ and approximately saturates above p± ~ 6GeV/c for y/s = 
200 GeV. Thus at high momenta the major fraction of 'hadrons' is of "leading" 
origin, i.e. pre-hadronic states containing quarks, antiquarks or diquarks from 
the primary struck nucleons, and may interact according to the assumptions 
stated above without delay with the (pre-) hadronic environment (cf. Fig. 2c). 

2.3 Initial state Cronin enhancement 

As known from the experimental studies of Refs. [46,47] an enhancement of 
the high transverse momentum particle cross section from proton-nucleus col- 
lisions - relative to scaled p+p collisions - is already observed at SPS and 
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ISR energies. This 'Cronin effect' is presently not fully understood in its de- 
tails, but probably related to an increase of the average transverse momentum 
squared (kj_) of the partons in the nuclear medium. One may speculate that 
such an enhancement of (k\) is due to induced initial semi-hard gluon radia- 
tion in the medium, which is not present in the vacuum due to the constraint 
of color neutrality. A related interpretation has been given by e.g. Kopeliovich 
et al. [17,35]. Since the microscopic mechanisms are beyond the scope of our 
present analysis, we do not want to comment this any further and simulate 
this effect in the transport approach by increasing the average (kj_) in the 
string fragmentation with the number of previous collisions A^ prev as 

(*1> = (kl) pp (l + aN pmv ) . (3) 



The parameter a ~ 0.25 — 0.4 in (3) is fixed in comparison to the experimental 
data for d+Au collisions [13,14] (see below). The assumption (3) is in line with 
an independent suggestion in Ref. [48]. 

In order to show the effect of enhancing the average (k\) in the fragmentation 
of the string we show in Fig. 4 the ratio of the charged hadron spectra from 
pp collisions for different (k\) to the spectra for the default value of {kj_) pp = 
0.36 GeV 2 . These ratios may be denoted as 'Cronin enhancement factors', 
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Fig. 4. The Cronin enhancement factor (see text) as a function of the transverse mo- 
mentum p± of charged hadrons for different average according to the Pythia 
calculations for N+N collisions at y/s = 200 GeV. 

which are a function of (k\) and the hadron species. We note, that these 
'Cronin enhancement factors' enter the dynamical transport calculations for 
each hadron type separately, whereas only the average enhancement factor is 
displayed for charged hadrons in Fig. 4. 
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The effect of such Cronin enhancement is obvious from Fig. 4: at low transverse 
momentum one observes a suppression of the p± spectrum with increasing 
(k\), which leads to a hardening of hadron spectra from nucleus-nucleus colli- 
sions relative to p+p reactions for transverse masses below about (mi-m ) < 
1 GeV. A maximum in the enhancement factor arises at p± ~ 2 ... 3 GeV/c 
followed by a smooth decrease towards unity at high transverse momenta, but 
a sizeable (constant) enhancement (above unity) remains up to the highest 
momenta considered in our calculations. 

2.4 d+Au collisions 

Whereas elementary p+p collisions give pseudorapidity distributiuons dN x /drj 
of produced hadrons x symmetric in 77 (Fig. 5), the recent measurements 
of dN x /dr) of primary charged hadrons in minimum bias d+Au collisions at 
x/s = 200 GeV from the PHOBOS collaboration [50] show a significant asym- 
metry in 77 (Fig. 5). This asymmetry in 77 is quite well reproduced by our 
calculations for minimum bias d+Au collisions (Fig. 5) as well as the approx- 
imate maximum in dN/drj for 77 w —2. This implies that the rescattering of 
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Fig. 5. The calculated pseudorapidity distributions dN/dr] of charged particles from 
inelastic p+p and d+Au reactions at y/s = 200 GeV (solid lines) in comparison with 
the data for p+p from the UA5 Collaboration [49] and d+Au collisions from the 
PHOBOS Collaboration [50]. [The experimental systematic errors for d+Au are on 
the level of 20%, whereas those in the normalization of our transport calculations 
are less than 10 %.] 

pre-hadrons from the deuteron projectile (77 > 0) with nucleons of the tar- 
get (77 < 0) is essential and that particles with large longitudinal momentum 
(77 > 2) are absorbed substantially in the Au-target ending up at pseudora- 
pidity ranges closer to the target pseudorapidity. 
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In order to sustain these findings, Fig. 5 shows also the comparison of our 
model calculations for p+p with data for p+p from UA5 [49] at the same 
invariant energy. The difference between the data and the calculations is on 
the 10% level and partly due to the fact, that UA5 measures inelastic events 
for p+p, whereas the calculations are for p+p collisions: We attribute the 
differences for pseudorapidity regions close to those of the initial baryons (an- 
tibaryon) to an excess of light quarks in the p+p reaction. 

Thus, in addition to the comparison in Fig. 1 for the transverse momentum 
distribution also the (soft) longitudinal hadron production is reasonably well 
described in our transport approach for elementary p+p as well as d+Au 
reactions. 

Since the Cronin enhancement factor essentially depends on (kj_), which in 
turn is controlled by the parameter a in eq. (3), we will determine its range 
by the data on d+Au reactions. 

A comparison of the calculated ratio 



is shown in Fig. 6 with the respective data for charged hadrons from the 
PHENIX and STAR collaborations [13,14] using a = 0.25 •••0.4 (hatched 
band). In eq. (4) (AT coll ) 8.5 denotes the number of inelastic nucleon-nucleon 
collisions per event (for minimum bias collisions), whereas 0"™ elas is the inelastic 
p+p cross section 1 . We find that our calculations give a rise in the ratio eq. (4) 
for small p± and a slight decrease for pj_ > 2.5GeV/c which can be traced 
back to inelastic interactions of the leading pre-hadrons with the nucleons in 
the target. The description of the data for RdA(p±) is not perfect and the data 
even indicate a slightly larger enhancement for p± = 4. . .6GeV/c, but the 
agreement is sufficient to proceed with Au+Au collisions at the same y/s. The 
important issue to remember in this context is that no dramatic absorption 
of high p± hadrons in d + Au collisions is found in the calculations as well as 
in the data. 

We note that the number of inelastic collisions (N co u) entering the ratio eq. (4) 
are directly available in the transport calculations; thus no model dependence 
enters the computation of this ratio. 

For completeness we show in Fig. 5 also the result for the ratio (4) from a 
calculation without any initial state Cronin enhancement (solid line). In this 

1 We remark with some caution, that the underlying binary scaling according to 
the Glauber model has recently been questioned by Kopeliovich [51] since inelastic 
shadowing corrections should reduce the total inelastic cross section by up to 20%. 



RdA(p_i) 



l/N^ d 2 N dA /dydp ± 




(4) 



(A^coii) /< clas d 2 a pp /dydp ± 
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Fig. 6. The suppression factor RdA for minimum bias d+Au collisions 
(•^/s = 200 GeV) at midrapidity for charged hadrons. Experimental data are from 
PHENIX [13] (open circles) and STAR [14] (filled squares). The hatched band de- 
notes our calculations for the 'Cronin' parameter a = 0.25 • • • 0.4 in (3) while the 
solid line results from transport calculations without employing any initial state 
Cronin enhancement (a = 0). 

limit the ratio drops below unity due to the interactions of the leading pre- 
hadrons with the residual target nucleons. 



2.5 Au+Au collisions 

As in case of the d+Au system we expect also an attenuation of high mo- 
mentum hadrons in Au+Au reactions due to inelastic interactions of leading 
pre-hadrons with the nucleons from the target (projectile). In addition, in- 
elastic interactions between the pre-hadrons among each other should occur 
as well as between leading pre-hadrons and formed secondary hadrons. Also in- 
teractions of formed secondary hadrons between each other might contribute 
to the final attenuation of high particles, if they are formed sufficiently 
early inside the hot/dense fireball [21]. We mention, that in the calculations 
to be shown below there are no interactions between leading pre-hadrons and 
explicit quark or antiquarks or gluons in the very early phase of the collision. 

Furthermore, in the HSD approach the formation of secondary hadrons is not 
only controlled by the formation time T/, but also by the energy density in the 
local rest frame, i.e., hadrons are not allowed to be formed if the energy density 
is above 1 GeV/fm 3 [28]. This energy density cut in HSD prevents hadrons 
to be formed in central nucleus-nucleus collisions at RHIC for a couple of 
fm/c when setting the clock by the initial impact of the two ions. We note 
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in passing, that a reduction of this energy density cut to 0.7 GeV/fm 3 (cf. 
Introduction) leads to a further delay of secondary hadron formation by ~ 1.2 
fm/c. Recall, however, that this energy cut does not apply for the pre-hadronic 
states defined above. 

We start with most central (5% centrality) Au+Au collisions at y/s = 200 GeV. 
In this case the nuclear modification factor is defined in accordance with eq. (4) 
as 

d /„ \ l/gXT ! d 2 N AA /dydp ± 

n A A{P±) (Arcou) /a ^ elas d 2 app/dydp± ■ W 



Fig. 7 shows a comparison of the calculations for eq. (5) with the data for 
charged hadrons from Ref. [10,11], where the hatched band again corresponds 
to the range in the parameter a in eq. (3) from 0.25 . . . 0.4. The calculations 
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Fig. 7. The suppression factor -Raa of charged hadrons at • • • 5 % central Au+Au 
collisions (y/s = 200 GeV) at midrapidity; Experimental data are from Refs. [10,11]. 
Please note, that the PHENIX data are for • • • 10 % centrality. The hatched band 
denotes our calculations for the 'Cronin' parameter a = 0.25 • • • 0.4 in (3). The solid 
line results from transport calculations without employing any initial state Cronin 
enhancement (a = 0). 

roughly reproduce the shape of the ratio -Raa(rl) but overestimate the ex- 
perimental ratio at higher p±. We emphasize, that the Cronin enhancement 
is visible at all momenta, but does not show up to be responsible for the 
peak structure in the enhancement around 2GeV/c. This becomes clear by 
comparing to the lower solid line, which results from transport calculations 
without employing any initial state Cronin enhancement. In this case the ratio 
-Raa(rl) is slightly low at p± ~ 2 ... 3 GeV/c in comparison to the data from 
the STAR collaboration, but rather well in line with the data from PHENIX. 
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For p± > 5GeV/c, however, it still underestimates the suppression seen by 
both collaborations. 



In order to understand this result we decompose the ratio Raa into contribu- 
tions from pions, kaons and p, u, K* vector mesons (before decays) without 
employing any initial state Cronin enhancement. The resulting comparison is 
shown in the left part of Fig. 8 for 5 % central Au+Au collisions. Here the 
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Fig. 8. The nuclear modification factor for central Au+Au collisions 
(i/s = 200 GeV) prior to any particle decays. The l.h.s. shows the total theoreti- 
cal ratios for the particle classes w, p, uj and K, K*. The right panel separates the 
theoretical contributions of the left panel into the 'leading' and the 'non-leading' 
particle contributions. 

direct pions show a reduced attenuation, the kaon reduction is slightly larger 
for lower p±, while the vector meson absorption is much stronger. Hadron for- 
mation time effects apparently play a substantial role in the few GeV/c region 
since heavier hadrons are formed earlier than light pions in the cms frame at 
fixed transverse momentum due to the lower Lorentz boost. 



This expectation is quantified in the right part of Fig. 8 where the decomposi- 
tion of the ratio Raa is additionally performed for leading and secondary /nonleading 
hadrons. First of all, the attenuation of the leading pions, kaons and vector 
mesons is roughly independent on p± and hadron type and gives Raa ~ 0.09. 
Thus the dominant attenuation seen experimentally should be addressed to 
the interactions of leading (unformed) pre-hadrons. This becomes even more 
apparent in accordance with all "back-of-the-envelope" estimates (including 
a constant rest frame formation time rf) when looking at the suppression fac- 
tors for the non- leading hadrons in Fig. 8 (right part). Here the direct pions 
practically do not show any attenuation in line with the expectation from 
the large formation times of pions with momenta of a couple of GeV/c. The 
situation changes for kaons and especially vector mesons due to their larger 
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mass. Now formed hadrons of a few GeV/c may interact with other formed 
hadrons and look 'thermalized', i.e. their final momenta show up in the soft 
part of the spectrum (below p± ~ 2 GeV/c). Since the expanding fireball has 
the same geometrical shape for all high p±_ hadrons, the attenuation in this 
"constant r/" picture is intimately correlated with the Lorentz 7-factor of the 
particles such that for high p± mesons above about 6 ... 8 GeV/c no absorp- 
tion of hadrons should occur via interactions of formed hadrons. On the other 
hand, the final attenuation of high p± hadrons is approximately independent 
on the actual formation time because according to Fig. 3 these particles are 
essentially leading pre-hadrons that rescatter anyhow. 

In order to demonstrate the effect of the formation time Tf we show in Fig. 9 
our results for the attenuation of secondary pions, kaons and vector mesons for 
Tf = 0.5fm/c. The latter value for 7/ is even favored by studies on antiproton 
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Fig. 9. The same as Fig. 8 (right panel) for a formation time Tf = 0.5fm/c. (Only 
the ratios of non-leading hadrons are displayed.) 

production and reabsorption in p+A reactions at AGS energies [45] as well 
as hadron attenuation in virtual photoproduction off nuclei [43]. In line with 
the discussion above we see from Fig. 9 - in comparison to the right part 
of Fig. 8 - that the suppression of kaons and vector mesons increases for a 
lower formation time and also extends to higher p±. However, as noted above, 
the total attenuation is essentially determined by the fraction of leading pre- 
hadrons at fixed p± (Fig. 3) and by their effective cross section with baryons 
or mesons. Formation time effects can only show up at lower p±. 

We note in passing that a decrease in the energy density cut from 1 GeV/fm 3 
to 0.7 GeV/fm 3 in the HSD calculations with respect to the formation of 
secondary hadrons leads only to a minor modification (enhancement) of the 
ratio Raa{p±) by about 10%, i.e. for p± > 6GeV/c the ratio increases from 
~ 0.4 to ~ 0.45. 
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We now turn to an alternative model for the leading pre-hadron cross section, 
which is discussed in the literature [52] (cf. also Ref. [21]) since the notion of a 
fractional constant cross section might be questionable [35] and alternative as- 
sumptions should be tested. 2 To this aim we have adopted a time-dependent 
cross section for leading pre-hadrons of the kind 

CleacKv^) = - 7 — Ohad(Vg) (6) 



for t— to <tf, where t denotes the actual production time and tf its formation 
time in the calculational frame. The full hadronic cross section is adopted for 
t > to- The numerical results of this assumption are shown for the ratio Raa in 
Fig. 10 in case of 5 % central Au+Au collisions. We find that the concept eq. (6) 
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Fig. 10. Same as fig. 7, but with a leading cross section according to eq. (6) for the 
perturbative high p± particles. 

is not in accordance with the experimental observation since the leading pre- 
hadron cross section is initially too low, i.e. when the fireball is very compact 
and dense. Furthermore, assuming eq. (6) in the nonperturbative transport 
calculations leads to dramatic consequences: In this case the pion, kaon and 
antibaryon rapidity distributions are severely underestimated due to the lack 
of inelastic reactions from leading pre-hadrons during the passage time of the 
heavy nuclei. 

We thus continue our analysis with the 'default' formation time Tf = 0.8fm/c 
and leading pre-hadron cross sections in line with the constituent quark model 



2 In the string fragmentation picture (combined with the constituent quark model) 
a considerable separation in space for the string ends at the starting time is un- 
avoidable. The only possible pointlike structure here is the hard scattered quark 
from a string end. 
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without attempting to 'fine-tune' any parameter. 
2.6 Dependence on centrality 

The centrality dependence of the ratio Raa is shown in Fig. 11 for 20 • • • 30 %, 
30 • ■ • 40 %, 40 • • • 60 % and 60 • • • 80 % centrality of Au+Au collisions at = 
200 GeV. Again the hatched bands correspond to our calculations with a 
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Fig. 11. Same as Fig. 7 for the centralities 20 • • • 30 %, 30 • • • 40 %, 40 • • • 60 % and 
60---80%. Calculational fluctuations are due to limited statistics. Experimental 
data are from Refs. [11,10]. 

Cronin parameter a in (3) ranging from 0.25 to 0.4 while the data stem from 
Refs. [11,10]. Note, that the uncertainty in the Cronin enhancement (width of 
the hatched band) decreases for more peripheral reactions, which is a direct 
consequence of the lower number of hard iViV collisions in (3). Our calcu- 
lations turn out to be in a better agreement with the data - or interpolate 
between the measurements from the STAR and PHENIX collaborations for 60 
to 80% centrality - than for the most central Au+Au collisions discussed in the 
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previous subsection. An analysis with respect to various contributions as in 
Fig. 7 shows, that with decreasing centrality the influence of formed hadrons 
essentially vanishes in the ratio Raa{p±) and the final state suppression is 
dominated by the interaction of leading pre-hadrons either 'immediately' with 
the baryons or other leading debries of the colliding nuclei or with the hadrons 
formed at later times. 

In order to obtain a measure for the centrality dependence of the suppres- 
sion we display in Fig. 12 the suppression factor (Raa) for charged hadrons 
integrated for p±_ > 4.5 GeV/c in comparison to the same quantity from the 
PHENIX collaboration [10] as a function of the number of participating nu- 
cleons iVpart- We mention that the determination of N paTt is model dependent 
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Fig. 12. The ratio (Raa) for charged hadrons for p± > 4.5GeV/c as a function of 
IVpart in comparison to the data from Ref. [10]. 

for very peripheral reactions, but rather save for mid-central and central re- 
actions. As already seen from Figs. 10 we describe the average suppression 
factors well for iVp art < 150 and underestimate the suppression for central 
Au+Au interactions. 

It is worth pointing out that the effective cross section of the leading pre- 
hadrons is constant during their propagation and thus the suppression is linear 
in the effective propagation length L. As shown in a geometrical model by 
Drees et al. [53] this linear absorption mechanism - when fixed in strength 
to the most central collisions - overestimates the amount of suppression for 
mid-central reactions. According to Ref. [53] the absorption geometry favors 
a suppression ~ L 2 in comparison to the data in Fig. 12. Though we do not 
find a strong argument in favor of any ~ L 2 suppression mechanism for all 
centrality classes we, nevertheless, seem to need some additional attenuation 
for central Au+Au reactions which could be attributed to strong interactions 
of the leading pre-hadrons in a 'colored' medium. 
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2. 7 Elliptic flow and Angular Correlations 



Since the leading pre-hadron cross section employed in our transport calcula- 
tion is much larger than related cross sections from pQCD, shadowing phe- 
nomena, that are related to the initial geometry in coordinate space for fixed 
impact parameter b, should show up in angular correlations like the elliptic 
flow of the high p± hadrons. In fact, a rather high and approximately con- 
stant elliptic flow has been seen experimentally for pj_ > 2GeV/c. It is still 
to be answered, whether this elliptic flow is a consequence of suppression 
due to strongly interacting matter or a survivor of correlations from the hard 
nucleon-nucleon interactions and following fragmentations. 

Gating on particles with momenta p± > 4 GeV/c in the transport calculation, 
only 1/3 of these final hadrons have suffered one or more interactions dur- 
ing their propagation to the vacuum, whereas the other 2/3 escape without 
any interaction. This observation implies, since more than 3/4 of the high p± 
hadrons are strongly absorbed, that the final high px hadrons seen experimen- 
tally essentially stem from pre-hadrons that originate from a diffuse 'surface 
region' of the expanding fireball. The latter pre-hadrons then evolve (or frag- 
ment) to the final hadrons dominantly in the vacuum and are accompanied by 
secondary hadrons, which - due to their large formation time - also hadronize 
in the vacuum. Thus, when gating on high p± hadrons (in the vacuum) the 
'near-side' correlations should be close to the 'near-side' correlations observed 
for jet fragmentation in the vacuum. According to our understanding the ab- 
sence of 'near-side' jet broadening does not necessarily indicate the presence 
of radiative quark energy loss in a QGP phase as advocated in Ref. [54]. 

We emphasize, that within our calculations the particles with large transverse 
momenta - that are treated perturbatively - so far have no intrinsic corre- 
lations as characteristic for back-to-back jet correlations. This point will be 
addressed in the near future within a separate study. 

The elliptic flow defined by 



provides a measure for the coupling of degrees of freedom in coordinate space 
with those in momentum space. In eq. (7) the y-direction denotes in-plane 
motion whereas the x-direction denotes out-of-plane propagation. At RHIC 
energies a strong in-plane elliptic flow is seen experimentally which for p± < 
1.5GeV/c is rather well described by hydrodynamical calculations [55]. How- 
ever, hydrodynamical calculations overestimate the approximately constant 
elliptic flow for higher transverse momenta. This observation proves that the 




(7) 
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medium created in Au+Au collisions is strongly interacting and the cross sec- 
tion needed to generate such a high elliptic flow is much larger than expected 
from standard kinetic pQCD [56]. 

In Fig. 13 we show the calculated elliptic flow (V2) for charged hadrons with 
transverse momentum p± = 2---6GeV/c in comparison to the data from 
Ref. [57]. One observes that our calculations underestimate the data by more 
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100 200 300 
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Fig. 13. The elliptic flow (i^) for charged hadrons for p±_ = 2---6GeV/c as a 
function of N part in comparison to the data from Ref. [57]. 

than a factor of 3; however, it is presently not clear if non-flow contributions 
of different origin contribute to the measured v 2 values. The elliptic flow de- 
scribed in our approach is essentially due to a shadowing of high p±_ hadrons 
in the out-of-plane direction, since there are no contributions from a collec- 
tive pressure in the perturbative dynamics of the high p± hadrons. As already 
addressed by Shuryak [58], such shadowing effects should be too weak to de- 
scribe the large v 2 values seen experimentally. The authors of Ref. [53] come 
to a very similar conclusion in their geometrical model (cf. also Ref. [54]). 

2.8 SPS energies 

In this subsection we briefly comment on results obtained in central collisions 
of Pb+Pb at -y/i = 17.3 GeV, i.e. at top SPS energies. We mention - without 
explicit representation - that the midrapidity p± spectra of hadrons from our 
Pythia calculations are approximately exponential in the transverse momen- 
tum showing no explicit power-law shape as for = 200 GeV in Fig. 1. This 
result stems in part from the kinematical cut in the maximum transverse mo- 
mentum, which is about 7GeV/c at ^/s = 17.3 GeV since two strings (jets) 
have to be formed. Furthermore, a direct comparison to related experimental 
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spectra is not possible because such data have not been measured. One might 
extrapolate the hadron p± spectra from higher energies assuming a x± scaling 
law to get a phenomenological parametrisation as in Ref. [59], however, our 
independent attempt in this direction showed systematic, non-constant errors 
(in important parts of the spectra) up to a factor of 2 when employing this 
ansatz. This uncertainty has to be kept in mind especially when comparing p± 
spectra from central Pb+Pb collisions to scaled 'extrapolated' spectra from 
p+p reactions at the same \fs. 

Since at yfs = 17.3 GeV the kinematical limits become essential for high p± 
hadrons, the ratio of 'leading particles' to 'all produced particles' (before par- 
ticle decays) in N+N collisions changes considerably in comparison to Fig. 3 
for y/s = 200 GeV. This information is displayed in Fig. 14 for different par- 
ticle classes as a function of transverse momentum as calculated within our 
Pythia description. A direct consequence of the results in Fig. 14 is that 
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Fig. 14. Same as Fig. 3 for y/s = 17.3 GeV. (Note the different momentum scales.) 
The thin dotted line around 0.1 shows the contribution of leading baryons with 2 
leading quarks. 

above ~ 2 ... 3 GeV/c now most of the mesons and the major fraction of nu- 
cleons are of leading origin and - within the dynamics specified above - may 
interact without delay. 

A further consequence of the steeper transverse momentum spectra at y/s = 
17.3 GeV relative to yfs = 200 GeV (Fig. 1) is that we get a larger Cronin 
enhancement for the SPS energy. This effect is due to a more drastic change in 
the spectra when enhancing (k\) in the string fragmentation; as a consequence 
the ratio between a calculation with enhanced (k\) to a calculation with the 
default value for p+p collisions becomes larger as seen from Fig. 15. This 
implies that the initial state Cronin enhancement is larger than at RHIC 
energies (cf. Fig. 4). 
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Fig. 15. Same as Fig. 4 for = 17.3 GeV. (Note the different momentum scales 
and also the different sizes.) 

Fig. 16 shows the result of our transport calculations for the ratio (5) in 5 % 
central Pb+Pb collisions at \fs = 17.3 GeV, where the hatched band again 
corresponds to the range in the parameter a in eq. (3) from 0.25- • -0.4. As 




P T [GeV] 

Fig. 16. The suppression factor -Raa of charged hadrons at • • • 5 % central Pb + Pb 
collisions (y/s = 17.3 GeV) at midrapidity. The hatched band corresponds to our cal- 
culations with the parameter a in (3) ranging from 0.25 to 0.4. The solid line results 
from transport calculations without employing any initial state Cronin enhancement 
(a = 0). 

in Fig. 6 the solid line shows a calculation without any initial state Cronin 
enhancement. We see from Fig. 16, that the uncertainty in the ratio (5) due 
to missing constraints on the initial state effects becomes large. The under- 
lying physics, however, is the same as at RHIC energies: the leading high p± 
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pre-hadrons interact with the baryons of the projectile/ tar get as well as with 
secondary formed hadrons. This leads to a net ratio of ~ 0.35 for p± >2 GeV/c 
when discarding any initial state enhancement (solid line in Fig. 16). However, 
this large attenuation is counterbalanced by a strong initial state Cronin en- 
hancement, which leads to an increase of the ratio (5) up to 2.5 ... 3 GeV/c 
followed by a moderate decrease. 

We have to mention that the WA98 collaboration has reported on the ratio 
Raa for neutral pions at top SPS energies [59] using extrapolations for the 
7T° spectra from p+p reactions taken from it (and 7r ± ) spectra at larger ener- 
gies. They quote a steep increase of -Raa with p±_ up to transverse momenta 
of 3.5GeV/c (although with large error bars), which is not described by our 
transport calculations. However, in view of the uncertainties in the normal- 
ization of the ratio (and the ratio itself!) we have to wait for measurements at 
RHIC for yfs w 20 GeV, where the transverse momentum spectra from p+p, 
d+Au and central Au+Au collisions can be measured by the same detector 
setup at the same y/s in order to clarify this issue from the experimental side. 



3 Summary 

Summarizing this study, we point out that (pre-) hadronic final state inter- 
actions, as included by default in present transport approaches like HSD, are 
able to approximately reproduce the high pj_ suppression effects observed in 
d+Au and Au+Au collisions at RHIC (y/s = 200 GeV). This finding is re- 
markable, since the same dynamics also describe the hadron formation and 
attenuation in deep-inelastic lepton scattering off nuclei at HERMES [43] 
quite well. Additionally, this also holds for antiproton production and atten- 
uation in proton-nucleus collisions at AGS energies [45]. 

We have, furthermore, found that interactions of formed hadrons after a for- 
mation time tf (which should be a matter of further debate) are not able 
to explain the attenuation observed experimentally for transverse momenta 
pj_> 6 GeV/c. However, the shape of the ratio Raa in transverse momentum 
p± reflects the presence of final state interactions of formed hadrons in the 
1...5GeV/c range [21]. 

In order to describe the experimental data below p± ~ 5 GeV/c, a reasonable 
description of the initial state "Cronin effect" has to be taken into account, 
which has been modeled in this work by increasing the average (k±) of the 
string fragmentation in subsequent hard collisions and has been fixed approx- 
imately by the d+Au measurements at RHIC. Such initial state effects are of 
pre-hadronic origin in a strongly interacting medium. Since they also show up 
in energetic proton-nucleus reactions [46,47], the initial state effects cannot 
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be considered as an indication for a phase transition (or crossover) of the high 
density matter produced in Au+Au collisions at relativistic energies. 

We have demonstrated, that the centrality dependence of the modification 
factor i?AA in Au+Au collisions at y/s = 200 GeV is well described for periph- 
eral and mid-central collisions on the basis of leading pre-hadron interactions. 
However, the attenuation in central Au+Au collisions is underestimated, which 
one may address to additional interactions of partons in a colored medium that 
have not been accounted for in our present transport studies. 

We finally note, that the elliptic flow i>2 for high transverse momentum parti- 
cles is underestimated by at least a factor of 3 (cf. Fig. 13). This demonstrates 
that geometrical shadowing in out-of-plane direction is not sufficient to get 
the right magnitude of the elliptic flow seen experimentally, which supports 
the independent studies in Refs. [58,53]. 

Since the high p± particles in the present work have no intrinsic back-to- 
back jet-like correlations, we postpone any further detailed statements to an 
upcoming investigation in the near future. Moreover, further experimental 
studies on the suppression of high momentum hadrons from d+Au and Au+Au 
collisions down to ^/s = 20 GeV will be necessary to separate initial state 
Cronin effects from final state attenuation and to disentangle to role of partons 
in a colored partonic medium from those of pre-hadrons in a color-neutral hot 
and dense fireball. 



The authors acknowledge valuable discussions with E. L. Bratkovskaya, D. 
d'Enterria, Th. Falter, B. Jacak, C. M. Ko, and Z. W. Lin. 
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